Symmetry breaking in frustrated systems: effective fluctuation spectrum due to coupling effects 
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We study the Langevin dynamics of a d-dimensional Ginzburg-Landau Hamiltonian with isotropic long range 
repulsive interactions. We show that, once the symmetry is broken, there is a coupling between the mean value 
of the local field and its fluctuations, generating an anisotropic effective fluctuation spectrum. This anisotropy 
has many interesting dynamical consequences. In the infinite time limit, static results are recovered which can be 
compared with the well known Brazovskii's transition to a state with modulated order. Our study reveals that the 
modulated solution appears continuously in d = 3 contrary to what is found in the usual approach neglecting 
coupling of modes, while for d < 3 transitions are still discontinuous. Analytical results for positional and 
orientational order parameters are also obtained and interpreted in the context of recent discussions. 
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Introduction.- According to a classic result due to Bra- 
zovskii m, systems in which the spectrum of fluctuations 
have a minimum in a shell in reciprocal space at a non zero 
wave vector, undergo a first order phase transition driven by 
fluctuations to a modulated state, in contrast to the second 
order transition predicted by mean field theory. This behav- 
ior, typically produced by the competition of isotropic interac- 
tions, have been successfully tested for two-dimensional mod- 
els in both theoretical |2] and numerical studies even 
in systems where interactions are not completely isotropic. 
Anyway, despite some efforts Isl j6t] there have not been such 
a clear confirmation for the character of the phase transition 
in other dimensions. The existence of a nematic phase 0- 
1^ have also opened the discussion about the nature of such 
a state with orientational but not translational order, that ap- 
pears above the Brazovskii transition. Since the strong degen- 
eracy of the fluctuation spectrum gives rise to the existence of 
many metastable structures, dynamical effects become very 
important |2, 6] and a great effort have been done [10l4l2ll 
in describing equilibrium properties by means of dynamical 
parameters. Following that way, in the present work we ad- 
dress some of the above points by solving the Langevin dy- 
namics of a standard model with competing isotropic inter- 
actions, breaking the symmetry by means of a small exter- 
nal field that is finally turned off. Taking the infinite time 
limit, we are able to study the steady state and recover the 
phase transition to modulated structures corresponding to that 
encountered in Brazovskii's equilibrium calculations flEltl- 
Moreover, we show that once the symmetry is broken, consid- 
ering non-homogeneous spatial contributions of the average 
order parameter on the fluctuations field, leads to the appear- 
ance of an effective fluctuation spectrum that is not isotropic 
anymore, modifying the usual interpretation of the low tem- 
perature phases. As a result, for = 3 we obtain a continuous 
transition to the modulated phase, contrary to what is gener- 
ally assumed, while the standard Brazovskii phenomenology 
is recovered when neglecting the coupling of modes. We find 



evidence of the existence of the nematic phase by evaluating 
analytical expressions for topological order parameters. This 
phase is interpreted in terms of the penetration depth of the 
fluctuations into the stripes structures, and is contrasted with 
a novel state also predicted, that consists in the presence of 
translational order without the orientational one. 

Model and Formalism.- We consider an effective Ginzburg- 
Landau Hamiltonian 
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where ?■ < 0, w > and J{x, x') = J{\x — x'\) represents a 
repulsive, isotropic, long range interaction. The external field 
is represented by h{x) and the parameter S measures the rela- 
tive strength between the attractive and repulsive parts of the 
Hamiltonian 1 10]. The Langevin dynamics takes the form; 
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d'^x' J{x,x')(j){x',t) + -^'q{x,t), (2) 



where F is set to 1 for simplicity, ri{x, t) represent a Gaus- 
sian noise that usually models the presence of temperature and 
h{x) — h cos(/co • x) represents a fictitious field that finally 
tends to zero and is introduced to break the symmetry. We 
may write the local field in terms of its fluctuations in the way 
4>{x,t) = ((/)) {x,t) + 'i>{x,t). Averaging over the thermal 
noise and initial conditions, the equation of motion (|2]l can be 
rewritten in terms of the average local field 
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where we have omitted for simplicity the dependence of mag- 
nitudes with {x, t) wherever it is possible. In the same way 
we can write from (|2]i and (O an equation for the fluctuations 
in the form 
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Here we have used the self-consistent Hartree approach, that 
is, the substitution of the term ( by 3 ( 5* in order to 
linearize the final equation. At this point the problem is to 
solve equations (|3]l and (|4]i for disordered initial conditions. It 
is worth to note that, from the linearity of equation (|4]i and the 
isotropy of the initial conditions, the quantity (4'^) must be 
independent of x. 

Stationary state.- Now, in the unimodal approximation, we 
assume {4){x,t oo)) — m cos(fco • x) where fco is the 
wave vector that minimizes the fluctuation spectrum. From 
equation (O it is direct to obtain the following result for m 
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where ^(fco) = 



(1/(5) J(fco) < 0, and the limit 



/i 0+ was made. Since (^t^) is an increasing function of 
temperature, equality (|5]l reveals the existence of a dynami- 
cal critical temperature Tc where real solutions cease to exist. 
The nature of this transition is not obvious because (^f^) must 
be obtained solving the whole set of equations, which may 
have no solution before m take the zero value. Now we have 
to solve equation (HJi in the infinite time limit; in the inverse 
space it can be rewritten 



(6) 
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where A{k) — A{k) — A{ko). In order to solve the equation 
system (|6]l we find the response function. Since there is a cou- 
pling between modes having a difference of momentum 2fco, 
all modes with a difference 2nko are actually coupled, this 
makes the exact solution to be a hard task. Instead, what we 
do is to consider the truncated system given by the two highest 
elements of the complete system's matrix. This approxima- 
tion captures the essential behavior of the highest eigenvalues 
of the system as functions of the momentum k. These eigen- 
values are the responsible of the long time behavior of the 
response function. We consider the highest eigenvalue that 
leads to a response function that does not diverge in the infi- 
nite time limit. Once we have the form of this response func- 
tion the calculation of {'i'^) can be done. The value of (^P^) 
depends on the topology of this eigenvalue as a function of k 
in the vicinity of fco- We use the form A{k) = ^(fc — fco)^ 
usually considered in similar contexts lilOil . In the infinite time 



limit we get: 
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where 7 = —A{kQ) — 2>u {^'^)- The expression (|7]) is the main 
result of this paper If we compare now with the usual kind of 
expressions obtained for (^^) we see that the system evolve 
like if it has an effective fluctuation spectrum of the form 



^eff(fc)cxA2(fci-fcof + ^(^|]fc,f^ 



(8) 



This determine a number of interesting features. The most 
obvious is the non-equivalence between spatial directions, this 
means that the properties of the system are not the same along 
the stripes and perpendicular to them. This particular form 
of the dispersion relation also has important consequences on 
the way in which transitions to the modulated phase occur for 
different dimensions. 

Mean square fluctuations.- In the context of the unimodal 
approximation 7 must be a small parameter, even more if we 
are near the transition. So, in order to study this transition we 
may expand integrals around 7 = 0. Defining the dimension- 
less groups: 



3u 
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where C{d) is an un-interesting constant depending on system 
dimension, we may obtain after some calculations 



(10) 



The solution of this equation gives us the value of (^I'^) for 
every temperature and, consequently, the critical temperature 
below which modulated structures appears. Contrary to the 
case d = 3 where the amplitude of the modulation goes to 
zero in a continuous way, for d = 1 and d = 2, before the 
annulation, this amplitude stops having a real solution. Solv- 
ing equation ( fTOl i for each dimension leads to the following 
critical values: 



d=l, 

d=2, 
d = 3, 




m (re) = ^ 

(V') (tc) = 0.777 
(V'')(rc)-1, 



(11) 



as can be seen in figure[T] There the behavior of ( V'^) is shown 
as a function of reduced temperature for every dimension. At 
this point we have obtained self-consistently the fluctuations 
for subcritical temperatures from a purely dynamical perspec- 
tive. We may now analyze the nature of the transition to the 
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FIG. 1: Temperature dependence of the amplitude of the principal 
mode m for different dimensions. The inset shows the corresponding 
values of the mean square fluctuations ('i/'^). 



modulated phase. From expression (|5]l it is possible to write 
the amplitude of the modulation m as a function of (-0^) and 
thus to obtain the value of m for every temperature in any di- 
mension. It is clear from figure[T]that, for d < 2, the transition 
to the high-temperature paramagnetic phase is discontinuous, 
contrary to d = 3, where m goes continuously to zero as r 
approaches to Tc- To better realize the consequences of the 
coupling between modes in our work let's consider the dy- 
namic calculation in the infinite time limit without the cou- 
pling terms in equation (|6]l. In this case the dynamical equa- 
tion is diagonal and modes with different momentum are not 
coupled. Following standard procedures, the self-consistent 
solution to the problem leads to the following equation: 



T 



(27r)° 



1 



7 + 
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where 7 — ~A{ko) — 3u (4'^) as before. What we have ob- 
tained is a new equation for the square mean fluctuations in 
which the fluctuation spectrum is spherically symmetric as in 
the standard Brazovskii-like approaches. The solution to this 
problem has been obtained many times in the literature |0, IJ. ] 



showing discontinuous transitions to the striped phase. Our 
new results suggest that the nature of the transition to the mod- 
ulated phase depends on dimensionality. 

Topological properties.- We now focus in the two- 
dimensional case in order to study the topological order of 
the modulated phase by means of two parameters commonly 
used in the literature with this purpose. ifisl [Till The first one 
is a positional order parameter defined by: 



M = 



d^x 4>{x) sin(A;o • x) 



(13) 



where S is the area of the system. As can be seen, this pa- 
rameter corresponds to a staggered normalized magnetization 
which is 1 only at zero temperature. Orientational order of the 



stripes can be quantified by means of the parameter 

Q^l [ d^x cos{2B{x)) 5{(t>{x))) / I I d^x 5{(t>{x)) 



(14) 

which again is defined in such a way that is 1 only in the per- 
fect stripes configuration. Here 5{4){x)) ~ 4){x + 5x) — (f>{x). 
In order to calculate parameters ( fT3] ) and ( fT4l) we use a gaus- 




FIG. 2: Behavior of A4 as a function of the reduced temperature for 
different values of c 

sian distribution for </> and its spatial derivatives in the evalua- 
tion of the corresponding averages. This assumption is consis- 
tent with the Hartree approximation and allows us to find ana- 
lytically closed expressions that, as far as we know, have never 
been obtained for topological order parameters. Since the final 
expressions are not particularly enlightning, we numerically 
evaluate them under some reasonable approximations. In this 
context we observe that the positional order parameter only 

depends on the dimensionless quantity c = ^S^°'' ■ In figure 
I— 1 

|2]we show typical behaviors of this positional parameter On 
the other hand, the orientational order parameter depends on 
c and also on the number A2. This dependence is shown in 
figure [3] As can be expected, the presence of a critical tem- 
perature for the two-dimensional system is reflected by means 
of jumps to zero of both, orientational and positional parame- 
ters at Tc (figures |2] and O. Our calculations show that these 
parameters decrease with increasing c, while the increase of 
A2 only causes an increasing of the orientational order This 
dependence of order parameters with A2 and c reveals a rich 
topological scenario for the two-dimensional system. In par- 
ticular, this picture supports the existence of the so-called ne- 
matic phase that has been the subject of some debate in the 
recent literature. 13] This nematic phase is characterized 
by a strong orientational order in the absence of positional or- 
der. Our results suggests that, by appropriately increasing c 
and A2, we may eventually obtain a region of temperatures 
close to Tc in which the orientational order prevails over the 
positional one. To test this one should have to consider a par- 
ticular model and its full fluctuation spectrum. On the other 
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FIG. 3: Behavior of Q as a function of the reduced temperature. The 
main figure corresponds to yl2 = 3 and several values of c. In the 
inset c = 5 and several values of A2 were taken. 

hand, the existence of a phase in which positional order pre- 
vails over orientational one is also suggested. From figure [3] 
one can see that increasing c the decay of the orientational 
parameter at low temperatures becomes more abrupt. This 
fact, together with the decrease of Q for increasing values of 
A2 (inset of figure 13) could be responsible for the appearance 
of such a phase. One can visualize this state by considering 
highly irregular stripe borders confined in such a small width 
that the staggered magnetization M is not too compromised. 
In that sense the lack of positional order in the nematic phase 
could be seen as the situation in which this width, that defines 
the penetration depth of the fluctuations, completely covers 
the width of the stripes. In this situation, thermal fluctua- 
tions are strong enough to break the stripes order generating a 
configuration with topological defects that, as have been dis- 
cussed in previous works, may give rise to a nematic phase. 
Finally, there is an absence of numerical evidence in literature 
supporting the existence of a phase with positional but not 
orientational order. This is probably due to the presence of a 
strong anisotropy, caused by a discretization mesh of the order 
of the modulation length in numerical works id, [Hi- In fact, 
this anisotropy effects favors the emergence of orientational 
order that makes the nematic phase easier to find. 

Conclusions.- Taking the infinite time limit, the steady state 
properties of a model of competing interactions have been 
studied from a purely dynamical calculation. Such a proce- 
dure makes use of the Hartree approximation over the fluctu- 



ations field and the isotropy of the system is explicitly bro- 
ken. This leads to an effective fluctuation spectrum that is not 
isotropic anymore, contrary to the one present in Brazovskii- 
like calculations. Within this formalism, it was possible to 
analytically obtain the fluctuations as functions of the tem- 
perature and thus analyze some interesting parameters in the 
modulated state. In particular, the first order character sug- 
gested in literature for this transition is not verified for the 
tree-dimensional case, in which a continuous transition takes 
place due to the coupling between modes. The present ap- 
proach extends our understanding about the nature of modu- 
lated phases in systems with competing interactions. We have 
also calculated analytically and discussed usual topological 
observables. Our restuls suggest that, for certain range of tem- 
peratures close to the transition, it may exist a nematic phase 
with purely orientational order, or even a phase in which the 
translational order prevails over the orientational one, depend- 
ing on some system dependent parameters. A more detailed 
paper including the discussion of dynamical observables will 
be published in short. 
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